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We show that, if graphene is subjected to the potential from an external superlattice, a band 
gap develops at the Dirac point provided the superlattice potential has broken inversion symmetry. 
As numerical example, we calculate the band structure of graphene in the presence of an external 
potential due to periodically patterned gates arranged in a triangular or a square graphene super- 
lattice (TGS or SGS) with broken inversion symmetry, and find that a band gap is created at the 
original and, in the case of a TGS, the "second generation" Dirac point. This gap, which extends 
throughout the superlattice Brillouin zone, can be controlled, in principle, by changing the external 
potential and the lattice constant of the superlattice. For a square superlattice of lattice constant 
10 nm, we have obtained a gap as large as 65 meV, for gate voltages no larger than 1.5 V. 

PACS numbers: 



I. INTRODUCTION 

Ever since the synthesis of high-quality graphene 
there has been tremendous interest in the properties of 
this single-layer form of carbon. Graphene has a honey- 
comb lattice structure, with two atoms per primitive cell 
and a hexagonal Brillouin zone (BZ). The Fermi energy 
Ep of homogeneous, neutral graphene lies at the so-called 
Dirac point, which occurs at high symmetry points K in 
the BZ. In fact, there are two inequivalent Dirac points 
Kq and Kq with two distinct valleys of excitations. Near 
the Dirac point, the density of states is linear in \E — Ep\ 
and the spectrum of quasiparticle states is well described 
by the Dirac equation for massless fermions. Partly as 
a result of this electronic structure [iJ-Q], graphene has 
many unusual electronic properties, such as a unique type 
of quantum Hall efFecH&jjBj , ballistic conduction by mass- 
less Dirac fermions size-dependent band gapQ, 
large magnetoresistance[J, H, and gate-tuneable op- 
tical transitions [Toj. 

A number of workers have investigated the possibility 
of building graphene electronic circuits without physi- 
cally cutting or etching the graphene monolayer. A nat- 
ural way to accomplish this is by subjecting graphene to 
an external potential with a suitable superlattice period- 
icity, e. g., by applying appropriate gate voltages. Such 
superlattices have proven to be extremely successful in 
controlling the electronic structure of more conventional 
semiconducting materials (see, e. g., Ref. [HI). In these 
materials, the presence of the additional periodic poten- 
tial gives rise to superlattice electronic band structure, 
which has extra band gaps at high symmetry points in 
the superlattice Brillouin zone (SBZ). 

There have been several predictions of electronic ef- 
fects in graphene due to an external superlattice poten- 
tial. For example, with one-dimensional (ID) and two- 
dimensional (2D) superlattices, the group velocity of the 
low-energy charge carriers is anisotropically renormalized 
[T^ . while a corrugated graphene sheet is expected to 



show charge inhomogeneity and localized states [Ts'l . Su- 
perlattices in graphene can be realized experimentally 
by using periodically patterned gates. Using an electron- 
beam, adatoms with a superlattice patterns of period- 
icity as small as 5 nm have been achieved on freestand- 
ing graphene membranes [HI . Superlattice patterns have 
also been observed for graphene on metal surfaces [isl - 
[iTj . In a recent theoretical study, a triangular graphene 
superlattice (TGS) was considered, and a new class of 
massless Dirac fermions was predicted to occur at the M 
point in the SBZfl^. 

In the studies so far, the external periodic potential 
giving rise to the superlattice has had spatial inversion 
symmetry. As a result of this symmetry, the degener- 
acy of the conduction and valence band at the original 
Dirac point (GDP) is preserved. In the present work, 
by contrast, we consider the potentials of an external 
TGS and a square graphene superlattice (SGS) without 
inversion symmetry. We show that, because of the ab- 
sence of inversion symmetry, an energy gap opens up at 
both the original and the new Dirac points in the TGS 
and at the Dirac point in the SGS. The magnitude of 
these gaps can, in principle, be controlled by modifying 
the externally applied voltages. Thus far, we have been 
able to achieve a band gap as large as 65 meV extending 
throughout the SBZ. Because of the controllability of the 
band gap, this type of graphene superlattice represents a 
system in which the band gap can be chosen with consid- 
erable freedom by modifying the applied external poten- 
tial. Such a system, and particularly the control at the 
ODP, might be quite useful in realizing future graphene 
electronic circuits. 



II. FORMALISM 

We first describe the formalism we use to calculate the 
superlattice band structure for a non-centrosymmetric 
superlattice potential, mostly following the approach of 
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Ref. 118|. If the periodicity of the superlattice is much 
larger than the intercarbon distance a ~ 1.42 A, which is 
the case considered here, the intervalley scattering (be- 
tween Ko and Kq) can be neglected [l3|. We limit our 
discussion to the spectrum near one of the two inequiva- 
lent Dirac points, which we denote Kg, in the presence of 
a periodic external potential V{x,y). In pure graphene, 

and 



we use a pseudospin basis ( ^ 



where ( ^ ) and 



are Bloch sums of tt orbitals with 



wave vector Kq on the sublattices A and B, respectively, 
and k is the wave vector relative to the Ko point. The 
single valley Hamiltonian for the quasiparticles of pure 
graphene is[20j 



(1) 



where vq is the (isotropic) group velocity at the Dirac 
point and the a's are Pauli matrices. The eigenstates 
and energy eigenvalues are 



<k(r) = 



V2 V 



and 



E°{k) = shvok, 



(2) 



(3) 



where s = ±1 is the band index and 0^ is the polar 
angle of the wavevector k. Eq. ([2]) indicates that the 
pseudospin vector is parallel and antiparallel to k in the 
upper (s = 1) and lower (s = —1) bands, respectively. 

We wish to consider the spectrum of elementary ex- 
citations near one of the Dirac points in the presence 
of V{x,y). The total Hamiltonian then takes the form 
H = Ho+H', where 



H' = IVix,y), 



(4) 



/ being the 2x2 identity matrix. We consider a periodic 
external potential having either a triangular Bravais lat- 
tice with basis vectors bi = bx, b2 = b(^^x+^y^, 

or a square Bravais lattice with bi = bx, b2 — by. 
We also assume that V{x, y) varies slowly on the scale 
of a lattice constant (6 ^ a). Then the band struc- 
ture of the elementary excitations near the Kq point can 
readily be obtained by diagonalizing Hq + H' in a plane 
wave spinor basis[llj. The basis states are of the form 

X(s)exp[i(k + G) • r], where x(s) = 75 ( 



or 
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and G is a reciprocal lattice vector cor- 
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responding to the periodic external potential V(x,y). If 
we include N plane waves, we must thus diagonalize a 
2A'' X 27V matrix to obtain the band structure of elemen- 
tary excitations non-perturbatively near the Kq point in 
the presence of the external periodic potential. We use 
N = 625 for both the triangular and the square graphene 



superlattice. In our actual calculations, we have consid- 
ered an external potential produced by a periodic array 
of circular regions, within each of which the potential is a 
constant. These are easy to treat by the above approach, 
because the Fourier transform of the constant potential is 
available analytically. To break the inversion symmetry, 
we have added a second, smaller circular region within 
each primitive cell, as illustrated in Figs. 1 and 2. If the 
second region is not equidistant from the two adjacent 
larger circular regions, inversion symmetry is broken. 



III. RESULTS 

We first review the case of a triangular array of large 
circles, with a constant external potential within each 
circle. Elsewhere in the superlattice, the external po- 
tential is zero. The resulting band structure has been 
computed in Ref. [ll|. In the presence of the superlat- 
tice, there are several Dirac points: the ODP, which is at 
r of the SBZ, and "next generation" Dirac points, at M 
and K of the SBZ. Of these new Dirac points, only the 
ODP (with energy set to eV) and the new Dirac point 
(NDP) at M with energy 0.196 eV) have vanishing den- 
sity of statesfl^. Near both of these points, the density 
of states goes linearly to zero with energy. However, the 
group velocity near the NDP is highly anisotropic. 

We now consider the superlattice band structure (SBS) 
for graphene subject to the non-centrosymmetric poten- 
tial shown in Fig. 1. We choose the lattice constant of 
the external potential to be 6 = 10 nm, and the radius 
ri of the large circles to be 2.5 nm, so that the filling 
fraction of large circles is 0.226725. We take the smaller 
circular gates to have radius r2 — 1 nm, and to be cen- 
tered at R + 0.45x, where R is a Bravais vector of the 
superlattice. The resulting band structure is shown in 
Fig. 3 for the case Vi = 0.5V, V2 = -1.65V, where Vi 
and V2 are the voltages on the larger and smaller circles. 
Once again = in the rest of the lattice. A substantial 
direct gap opens up at the ODP of magnitude Ar = 58.2 
meV for these parameters (between the bands above and 
below the ODP). Although not obvious from Fig. 3, there 
is also a gap A^/ which opens up at the NDP, of mag- 
nitude ~ 1 meV. For these parameters, we have checked 
that there is, in fact, a full band gap at ODP extending 
through the entire SBZ. The small gap also at M does 
not extend throughout the SBZ, as is evident from the 
Figure (the gap at M does extend throughout the SBZ 
for V2 positive). The inset of Fig. 3 shows the SBZ, with 
the symmetry points F, K , and M indicated. 

In Fig. 4, we show the corresponding SBS for graphene 
subjected to the periodic potential with square symmetry 
shown in Fig. 2. Once again, we choose the superlattice 
constant 6 = 10 nm, ri — 2.5 nm, and r2 — 1 nm, so 
that the filling fraction of circles in this case is 0.227765. 
We choose Vi = 0.5 V and V2 = —1.4 V, and each small 
circle is located at 0.375i relative to the center of each 
large circle. There is a gap of magnitude 64.4 meV, at 
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the ODP at T which extends throughout the SBZ. 

In Fig. 5, we show show how the direct gaps at the 
Dirac point depend on the voltage V2 on the smaher cir- 
cles in the triangular and square superlattices. For the 
square superlattice, we take & = 10 nm, ri/b = 0.25 
r2/b — 0.125, Vi at 0.5 V, and the small circles are lo- 
cated at 0.375a;. For the triangular superlattice, & = 10 
nm, n/b = 0.25, r2/b = 0.1, Vi = 0.5 V, and the small 
circles are located at OAx. In both cases, these gaps show 
clear maxima as functions of V2 ■ We have carried out sim- 
ilar calculations using positive values of V2 in both cases 
and also obtain nonzero gaps extending throughout the 
SBZ. In the case of positive V2, they are generally smaller 
than for negative voltages. 

In Fig. 6, we show the dependence of the gaps at the 
Dirac point on the radius r2 in both the triangular and 
square superlattices. For the square superlattice, we take 
& = 10 nm, Vi = 0.5 V, V2 = -1.4 V, n/b = 0.25 and 
the small circles are located at 0.3752;. For the triangular 
case, 6 = 10 nm, Vi = 0.5 V, V2 = -1.65 V, n/b = 0.25, 
and we locate the small circular contacts at OAx. The 
direct band gaps can reach as large as ~ 65 meV for the 
square superlattice. 

Because of the special form of the Hamiltonian of 
eqs. (1) and (4), the superlattice band structure, and 
hence the band gap at the ODP, exhibits a simple scal- 
ing relation: under the transformation 6 — A6 and 
V{x,y) — V{Xx, Xy)/ X, where A is a dimensionless scal- 
ing factor, the eigenvalues of H satisfy E{kx,ky) — 
E{kx/X,ky/X)/X. Thus, in particular, Ar is multiplied 
by 1/A. For example, if b is halved (and all other lengths 
are also halved), and the potentials on all the contacts 
are doubled, then the band gap is also doubled. We have 
confirmed numerically that this scaling relation is satis- 
fied for our bands. 

Besides these calculations, we have tried some other ge- 
ometries in an effort to obtain the largest possible band 
gap. For example, we have calculated the SBS for trian- 
gular and square superlattices with two off-center circular 
contacts, one along each edge of the primitive cell. This 
geometry does not lead to band gaps larger than we have 
obtained with one off-center circle. 



IV. DISCUSSION 

The present results show that the band gap of graphene 
at the Dirac point can be manipulated by subjecting it 
to a superlattice potential which lacks inversion symme- 
try. The magnitude of the band gap can be controlled 
by modifying the contact potentials. Thus, if such an ar- 
rangement of contacts can be created, the resulting ma- 
terial has a band gap which can be controlled via the 
applied voltages without modifying the structure. 

There are several obvious challenges before this scheme 
could be used in practice. First, the band gaps are rather 
small (of order 65 meV). These can be made larger by 
reducing b, while simultaneously increasing Vi and V2 



as implied by the above scaling relation, but this could 
be challenging experimentally. The gaps can also be in- 
creased by increasing the voltages at fixed geometry. Per- 
haps the most promising scheme might be to increase b 
while simultaneously increasing the bias voltages. The 
larger b would be easier to achieve experimentally, while 
the large voltage offsets would increase the gap. These 
ideas certainly do not exhaust the possibilities offered by 
periodic arrangements of contacts. Any 2D superlattice 
of gates lacking inversion symmetry would lead to a band 
gap at the ODP, and some which we have not tried may 
lead to a larger band gap than those we have found to 
date. A ID superlattice potential lacking inversion sym- 
metry will not suffice to produce a complete band gap. 
But such a ID potential, when combined with a suitable 
non-time-reversal-invariant perturbation such as a mag- 
netic field, might also lead to a complete band gap. 

There are, of course, many other ways to create a band 
gap in graphene besides the method described here. For 
example, one could create a non-centrosymmetric lattice 
of nanoscale holes in graphene. It would be difficult to 
use this method, however, to create a tunable band gap. 
In graphene nanoribbons with armchair edges, the band 
structure becomes insulating if the width of the sam- 
ple, in units of the lattice constant a, is not of the form 
3M + 1, with M an integer [2l| - |25| . But creating such a 
ribbon requires cutting graphene samples with very high 
precision, and again, once the ribbon has been created, 
the band gap cannot be easily tuned. Other possibilities 



are to induce a 



Ja' 



ip by using a substrate which lacks inver- 
sion symmetry [26 1, or to use graphene antidot lattices [27j 
Again in this case, the gap cannot easily be controlled be- 
cause the atomic spacing of the substrate is fixed. Our 
proposal of a TGS (or other 2D superlattice) with a non- 
centrosymmetric superlattice potential is more efficient in 
opening up a gap at the ODP. In principle, moreover, this 
gap can easily be larger than the thermal energy at room 
temperature for modest values of potential. Thus, this 
method may be a viable approach to creating graphene 
with a readily tunable band gap. 

In summary, we have calculated the superlattice 
band structure for graphene subjected to a non- 
centrosymmetric superlattice of contacts, on each of 
which the voltage is held constant. We find that the 
superlattice band structure exhibits a band gap extend- 
ing throughout the SBZ, which can thus be tuned by 
external voltages. For some choices of the voltages and 
superlattice constant, the band gap can be as large as 
65 meV, significantly larger than the room temperature 
thermal energy. Thus, this arrangement might possibly 
be of use in future electronic or electromagnetic devices. 
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FIG. 1: (Color online.) Schematic of graphene plane sub- 
jected to a periodic external potential V{x, y) consisting of 
a periodic arrangement of circular contacts, within each of 
which V is a constant. The large circles are arranged on a 
triangular lattice of edge h. The small circles are not mid- 
way between the large circles; for the case shown, they are 
centered at QAbx relative to the centers of the large circles. 
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FIG. 2: (Color online.) Same as Fig. 1, except that the large 
and small circular contacts are arranged on a square lattice. 
Each small circle is now centered at 0.3756a; relative to the 
center of the nearest large circle. 
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FIG. 3: (Color online) Superlattice band structure (SBS) for 
graphene subjected to an external potential produced by a 
triangular superlattice of contacts, as shown in Fig. 1. The 
calculation is carried out using 6 = 10 nm, Vi = 0.5 V, V2 ~ 
— 1.65 V, ri/b = 0.25, r2/b = 0.1, and the small circles are 
located at 0.4a;. The F point of the superlattice corresponds 
to one of the K point of the original graphene lattice. If only 
the large circular contacts were present, the SBS would have 
Dirac points at both F and M of the SBS Brillouin zone, at 
both of which the density of states would go linearly to zero. 
Inset: first Brillouin zone of a triangular lattice showing the 
three high-symmetry points F, M, and K. 
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FIG. 4: (Color online) SBS for graphene subjected to an ex- 
ternal potential produced by a square superlattice of contacts, 
as in Fig. 2. The F point of the superlattice corresponds to one 
of the K point of the original graphene lattice. In this case, 
we use 6 = 10 nm, Vi = 0.5 V, V2 = -1.4 V, ri/fe = 0.25, 
r2/b = 0.125 and the small circles are located at 0.375a;. In- 
set: first Brillouin zone of a square lattice showing the three 
high-symmetry points V, X, and M. 
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FIG. 5: (Color online) Plot of the gap at the Dirac point F as 
a function of the voltage V2 (in V) on the smaller circles. Full 
curve: square lattice; dashed curve: triangular lattice. For 
the square lattice, we use & = 10 nm, Vi = 0.5 V, ri/b = 0.25, 
r2/6 = 0.125 and small circles are located at 0.375i;. For 
the triangular lattice, 6 = 10 nm, Vi = 0.5 V, ri/b = 0.25, 
r2/b — 0.1 and the small circles are located at OAx. 
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FIG. 6: (Color Online) Plot of the gap at the Dirac point 
as a function r-^jh. Full curve: square lattice; dashed curve: 
triangular lattice. For the square lattice 6 = 10 nm, Vi = 0.5 
V, Vi — —1.4 V, ri/fe = 0.25, and small circles are located 
at 0.375a;. For the triangular lattice, b = 10 nm, Vi = 0.5 V, 
V2 ~ —1.65 V, ri/b = 0.25 and small circles are located at 
OAx. 



